#Get Data
polys_sf <- readRDS('./clean_data/pixel_values/polys_sf.rds') %>% st_make_valid() %>% as_tibble()
val_pred <- readRDS("./clean_data/pixel_values/val_pred.rds") %>%
left_join(polys_sf) %>%
dplyr::select(c(label, .pred_class, poly_num, area))
## Joining, by = c("label", "quadrat", "poly_num")
#scaleFUN <- function(x) sprintf("%.2f", x)
val_pred %>% group_by(label) %>% summarize(area = sum(area), n = n()) -> area_summary
ggplot(area_summary, label = n) +
geom_col(mapping = aes(x = label, y = area, fill= label)) +
# geom_text(nudge_y = 100000, size = 3) +
scale_fill_manual(values=c("#b7a386", "#5a6233", "#5e4a31"))# +
#scale_y_continuous(labels = scaleFUN)
val_pred %>% filter(label == 'live vegetation') -> live_vegetation_data
val_pred %>% filter(label == 'dead vegetation') -> dead_vegetation_data
val_pred %>% filter(label == 'sand') -> sand_data